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Abstract. We present a quantitative study of massive stars in the High Excitation Blob N81, a compact star 
forming region in the SMC. The stellar content was resolved by HST and STIS was used to obtain medium 
resolution spectra. The qualitative analysis of the stellar properties presented in Heydari-Malayeri et al. {2002a I 
is extended using non-LTE spherically extended atmosphere models including line-blanketing computed with 
the code CMFGEN (Hillier & Miller 1998), and the wind properties are investigated. The main results are the 
following: 

The SMC-N81 components are young (~ 0-4 Myrs) O stars with effective temperatures compatible with medium to late 
subtypes and with luminosities lower than average Galactic O dwarfs, rendering them possible ZAMS candidates. 
The winds are extremely weak: with values of the order 10 -8 /10 -9 Mq yr _1 the mass loss rates are lower than observed 
so far for Galactic dwarfs. Only the recent study of SMC stars by Bouret et al. (2003) show the same trend. The modified 
wind momenta (M Voo are also 1 to 2 orders of magnitude lower than observed for Galactic stars. Both the mass loss 
rates and the modified wind momenta are lower than the predictions of the most recent hydrodynamical models. 

The accuracy of the UV based mass loss rate determination, relying in particular on the predicted ionisation 
fractions, are carefully examined. We find that M could be underestimated by a factor of up to 10. Even in this 
unlikely case, the above conclusions remain valid qualitatively. 

The reasons for such weak winds are investigated with special emphasis on the modified wind momenta: 

There may be a break-down of the wind momentum - luminosity relation (WLR) for dwarf stars at low luminosity (log 
L/L0 Si 5.5). However, reasons for such a breakdown remain unknown. 

The slope of the WLR may be steeper at low metallicity. This is predicted by the radiation driven wind theory, but the 
current hydrodynamical simulations do not show any change of the slope at SMC metallicity. Moreover, there are indications 
that some Galactic objects have wind momenta similar to those of the SMC stars. 

Decoupling may take place in the atmosphere of the SMC-N81 stars, leading to multicomponent winds. However, various 
tests indicate that this is not likely to be the case. 

The origin of the weakness of the wind observed in the SMC-N81 stars remains unknown. We suggest that this 
weakness may be linked with the youth of these stars and represents possibly the onset of stellar winds in recently 
formed massive stars. 
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1. Introduction 

Massive stars play key roles in various astrophysical con- 
texts all along their evolution: they ionise ultra-compact 
HII regions while still embedded in their parental molec- 
ular cloud; they create ionised cavities and shape the sur- 
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rounding interstellar medium during the main fraction of 
their lifetime; they experience strong episodes of mass loss 
when they become Luminous Blue Variables and Wolf- 
Rayet stars, revealing their core and enriching the ISM in 
products of H and He burning; they end their life as super- 
novae, producing the heavy elements and releasing large 
amounts of mechanical energy. During all these phases, 
massive stars lose mass through winds driven by radiation 
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pressure on metallic lines. This affects not only their evo- 
lution (e.g. Chiosi & Maeder [1986(1 but also the surround- 
ing interstellar medium in which the release of mechanical 
energy can trigger instabilities leading to the collapse of 
molecular clouds and to star formation. Moreover, bubbles 
and superbubbles observed on galactic scales are powered 
by such mass ejections. Hence, various astrophysical fields 
require the knowledge of quantitative wind properties of 
massive stars. 

Several studies have been carried out in the last two 
decades to determine these properties. At solar metallicity, 
the observational determinations (e.g. Howarth & Prinja 
115551 Puis et al. 1X5551 Herrero et al. l3tjUU|> are on av- 
erage in good agreement with the most recent hydrody- 
namical predictions based on the radiation driven wind 
theory (Vink et al. I2000|l . both in terms of mass loss rate 
and of the modified wind momentum - luminosity rela- 
tion (WLR, e.g. Puis et al. [1996) which quantifies the 
strength of the wind. At non solar metallicities, we expect 
the wind properties to vary with Z due to the modified ra- 
diative acceleration through metallic lines. In particular, 
the mass loss rate should be proportional to Z r (Abbott 
H5521 Puis et al. EUODj l and the WLR should be shifted 
towards lower values and should have a steeper slope. The 
most recent theoretical results predict r ~ 0.8 (Vink et 
al. 12001(1 but no change in the slope of the WLR, at least 
for Z > 1(T 3 Z Q (Hoffmann et al. [2TJU21 KudritzkiEHS- 
Observational studies indicate a reduction of the mass loss 
rate and of the terminal velocity in the Magellanic Clouds, 
but given the small number of objects studied so far, the 
behaviour of the WLR at low metallicity is still poorly un- 
derstood. Several groups are currently analysing stars in 
sub solar (Crowther et al.[2002 Hillier et al. ESHEH Bouret 
et al. 12003(1 and super solar (Najarro et al. 119971 , Figer 
et al. 12002(1 regions for a better understanding of wind 
properties in different environments. The present work on 
SMC-N81 stars takes part in this effort. 

The SMC-N81 region belongs to the class of the "High 
Excitation Blobs" (HEB) first introduced by Heydari- 
Malayeri & Testor (1982). These blobs are compact re- 
gions of star formation in the Magellanic Clouds (see 
Heydari-Malayeri 2002b for a complete review). They have 
a typical radius of a few pc and display the features of 
star forming regions: HII cavities, turbulent structures, 
ionisation fronts and shocks. Recent HST observations 
(Heydari-Malayeri et al. [D299) have revealed for the first 
time its stellar content, 4 of the brightest stars being 
grouped in the central 2 " wide region. Subsequently spec- 
tra of the main exciting stars have been obtained with 
STIS onboard HST. The qualitative analysis of these spec- 
tra, presented in Heydari-Malayeri et al. (2002a, here- 
after paper I), have already revealed interesting proper- 
ties. First, the stars have been identified as mid O dwarfs 
with surprisingly low luminosities compared to "classi- 
cal" dwarfs. Second, the UV spectra have shown signa- 
tures of very weak winds, even weaker than those usu- 
ally observed in the SMC. These characteristics have lead 
Heydari-Malayeri et al. il2002a ■ to propose that the SMC- 



N81 stars could belong to the class of Vz stars which 
are massive stars thought to lie very close to the ZAMS 
(Walborn & Parker 11552)1 . 

As such the properties of these stars, showing unusu- 
ally weak winds compared to other SMC O stars, seem 
already quite interesting. Furthermore the association of 
these objects with a compact star forming region, pre- 
sumably indicative of a very young age, allows one also 
to obtain unique constraints on properties of very young 
massive stars shortly after their birth. In fact such ob- 
servations appear crucial for a better understanding of 
the earliest evolutionary phases of massive stars and to 
constrain their formation process which is still under de- 
bate (they may form by accretion on a protostellar core - 
Norberg & Maeder 125551 Behrend & Maeder I255T1 - or by 
collisions between low mass components in dense stellar 
clusters - Bonnell et al. 11998(1 . 

With such objectives in mind we have carried out a 
quantitative study of the UV spectra of the SMC-N81 
stars. First results have been presented in Martins et al. 
(2002a). In fact we are able to determine upper limits on 
the mass loss rates of four O stars in this region, which 
turn out to be surprisingly low (typically M a few 10 -9 
Mq yr _1 ) compared to predictions of the radiation driven 
wind theory, even when taking metallicity effects into ac- 
count. Although no precise physical explanation is found 
for this behaviour we strongly suggest that this behaviour 
is related to the very youth of these massive stars. 

The remainder of the paper is structured as follows. 
Section |2J briefly summarises the observations and data 
reduction. Section describes the main ingredients of the 
modeling. In Sect. 01 we explain how interstellar lines are 
taken into account. The main results are given in Sect. El 
and discussed in Sect. (nebular and stellar properties) 
and (wind properties). Finally, Sect. |5J summarises the 
main results. 



2. Observations and data reduction 

Ten SMC N81 stars were observed with STIS onboard 
HST on 28 and 31 October 1999 (General Observer 
Program No 8246, PI M. Heydari-Malayeri). The spec- 
tra in the wavelength range 1120 - 1715 A where ob- 
tained through the G140 L grating on the Multi-Anode 
MicroChannel Array (MAMA) detector. The 52" x 0.2" en- 
trance slit was used. The effective resolution was 0.6 A per 
pixel of 25 /im which correspond to a dispersion of 24 A 
mm -1 or a resolution of 1.2 A. The exposure times were 
chosen to equalize the S /N ratios which are of the order 20 
for the 4 stars studied here. The other ones have lower S/N 
ratio which precludes any quantitative analysis. Optical 
spectra, also obtained with STIS for 6 objects, were of 
insufficient quality (too low spectral resolution) and are 
therefore not used here. 

Details concerning the data reduction can be found in 
paper I. 
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3. CMFGEN: basic concepts 

The modeling of realistic atmospheres of massive stars 
requires the inclusion of three main ingredients: 1) due 
to the high luminosity of these stars, radiative processes 
are dominant and a non-LTE treatment must therefore be 
done; 2) mass loss creates an atmosphere which can extend 
up to several tens of stellar radii which renders the use of 
spherical geometry unavoidable; 3) the inclusion of metals 
(mostly iron) is fundamental to reproduce realistic atmo- 
spheric structures and emergent spectra (line-blanketing 
effect). 

We are now in an area where powerful tools including 
most of the above ingredients with progressively fewer ap- 
proximations are becoming available. Examples are the 
codes TLUSTY (Hub eny & Lanz 113951 . WM-BASIC 
(Pauldrach et al. lMfTjl and FAST- WIND (Santolaya-Rey 
et al. 11997}! . For our study, we have chosen to use the pro- 
gram CMFGEN (Hillier & Miller HUE)) now widely used 
for spectroscopic studies of massive stars. The main ingre- 
dients taken into account in CMFGEN are the following: 

— non-LTE approach: the whole set of statistical equi- 
librium equations and radiative transfer equations is 
solved to yield the level populations and the radiative 
held. 

— wind extension: all equations are written in spherical 
geometry with the assumption of spherical symmetry 
and include all velocity terms due to the expanding 
atmosphere. 

— line-blanketing: metals are included in the statistical 
equilibrium equations so that accurate level occupa- 
tion numbers can be derived. A super-level approxi- 
mation consisting in gathering level of close energy in 
a unique super-level is used to reduce the computa- 
tional coast. 

— radiative equilibrium: the temperature structure in the 
atmosphere is set by the condition of radiative equilib- 
rium. Other heating/cooling sources are optional and 
can be included (adiabatic cooling, X-rays). 

— hydrodynamical structure: at present, CMFGEN does 
not compute the velocity and density structure so that 
they have to be taken as input data or have to be 
parameteriscd. In our approach, they are computed 
by other atmosphere codes. We use either the ISA- 
WIND code which uses the hydrodynamical equations 
together with a grey LTE temperature to yield the den- 
sity and velocity in the atmosphere (see also Martins 
et al. 12002b}) . or TLUSTY which, thanks to a more 
accurate description of the pressure terms, gives a 
good description of the photospheric structure which 
is connected to a wind velocity structure represented 
by a classical law (v = Voo(l — ^r) 13 )- Tests made 
with stellar and wind parameters typical of our stars 
have shown that both methods give similar results for 
dwarfs in terms of emergent UV spectra. 

Several input parameters have to be specified, the main 
ones being: 



Table 1. Ions included in the atmosphere models. Number 
in parenthesis indicate species which, for given T e g, are 
trace ions and then are not taken into account. 



Element 


Ionisation state 


H 


I II 












He 


I II 


III 










C 


(II) 


III 


IV 


V 






N 


(II) 


III 


IV 


V 


VI 







(II) 


III 


IV 


V 


VI 


(VII) 


Si 


(II) 


(III) 


IV 


V 






s 




(III) 


IV 


V 


VI 


(VII) 


Fe 




(III) 


IV 


V 


VI 


VII (VIII) 



— stellar parameters: luminosity (L), radius (R), mass 
(M). 

— wind parameters: mass loss rate (M , terminal velocity 
Woo)- Note that the slope of the velocity field in the 
wind (the so-called /? parameter) is chosen when the 
velocity structure is computed so it is also an input pa- 
rameter (with default value 0.8). CMFGEN also gives 
the possibility to include clumping. This is done by 
the inclusion of a volume filling factor / of the form 
/ = foe + (1 - /oo)e~^ where is the value of / 
in the outer atmosphere and Vi n it the velocity at which 
clumping appears (30 km s _1 by default in our com- 
putations). 

— abundances/elements: in most of our model, the met- 
alicity have been chosen to be 1/8 solar where solar 
refers to the values by Grevesse & Sauval l}19981 . This 
metallicity is thought to be typical of stars in the SMC 
and in N81 in particular (Venn H^l Hill lT§551 Vermeij 
et al. 12002(1 . although the exact value remains uncer- 
tain. The ions included in our computations are given 
in Tabled 

— turbulent velocity: a value of 20 km s _1 was chosen 
for the calculation of the populations and of the tem- 
perature structure. Martins et al. (}2002b|) have shown 
that reasonable changes of v tUT b na d few effects on this 
part of the calculation. For the formal solution of the 
radiative transfer equation leading to the final detailed 
emergent spectrum, the value of Vt U rb has been deter- 
mined to give the best fit as shown in Sect. 15.51 and is 
found to be of the order 5 km s -1 . 



4. Interstellar lines 

The determination of mass loss rates relies on the fit of 
emission or P-Cygni lines (see Sect. 15. 2f) . The synthetic 
profiles are quite sensitive to the value of M, so that a re- 
liable estimate of this quantity requires the best possible 
knowledge of the stellar and wind spectrum. Consequently, 
the contaminating interstellar (IS) lines must be identi- 
fied, which is easy in high resolution spectra where these IS 
lines appear as narrow absorptions in the P-Cygni profiles. 
However, in our medium resolution observations the IS 
components are diluted in the stellar features so that the 
exact stellar + wind profiles are uncertain. Of course, this 
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Table 2. Determination of CIV column densities 



Component log(iV(C IV)) Reference 



Galactic 


14.4/14.51^ 


Mallouris et al. (2001) 




14.4 


Fitzpatrick & Savage (1983) 




14.35lH 6 6 


Sembach & Savage (1992) 




14.28l° j° 


SMC stars, this work 


SMC 


14.4/14.51°:!° 


Mallouris et al. (2001) 




> 14.5 


Fitzpatrick & Savage (1983) 




14 4Q+0-20 


SMC stars, this work 



depends on the line: as N V is a trace ion in the interstel- 
lar medium, stellar N v A1238,1242 is weakly affected by 
this problem. But this is not the case of C IV AA1548,1551 
which is all the more contaminated given that the wind 
feature is weak in the observed N81 spectra. It is there- 
fore crucial to estimate the contribution of the interstellar 
absorption to derive reliable mass loss rates. 

The interstellar absorption originates both in the 
Galaxy and the local SMC environment. To estimate the 
interstellar C iv column densities, we have proceeded as 
follows: 

- First, we have used high resolution HST-STIS spec- 
tra of 8 SMC stars (AV 69, AV 75, AV 80, AV 327, 
NGC346 355, NGC346 368, NGC346 113, NGC346 12), 
for which UV spectra have been obtained from the HST 
archive, to determine the C iv column densities in the 
direction of the NGC 346 region and the southwest part 
of the SMC. For that purpose, the Galactic and SMC in- 
terstellar profiles have been fitted with a gaussian pro- 
file (with a shift of 140 km s^ 1 - the receding velocity - 
for the SMC component). Two parameters are needed to 
achieve such a fit: the column density and the FWHM of 
the gaussian profile (for which a typical value of 20 km s _1 
was chosen). The mean column densities derived from this 
study are the following: log(V(C IV)) = 14.28T.£^ for the 
Galactic component, and log(JV(C IV)) = 14.43±g;^ for 
the SMC component. The higher dispersion in the case of 
the SMC component is probably due to the fact that we 
are looking at different parts of the SMC, whereas only 
one line of sight is used for the Galaxy. 

- Second, we have taken various determination of the 
C IV column density from the literature. Several methods 
were used (curve of growth, line fitting, apparent optical 
depth) and give consistent results. 

Table |3 summarises the various column density esti- 
mates. For the Galactic case, the results from our deter- 
mination are consistent with more accurate determina- 
tions (within the errors bars). For subsequent analysis, 
we adopt log(iV(C IV) = 14.4 as a representative value 
for the Galactic absorption. For the SMC, the results are 
also in good agreement. As the redenning of the SMC-N81 
stars is similar to that of the other 8 SMC stars used for 
this study (the values of E{B — V) are of the order 0.12 for 
all stars) , the properties of the interstellar matter on the 
different lines of sight towards the SMC sampled here must 
be the same. This absence of local extinction was noted by 
Heydari-Malayeri et al. {1988 ) in the first detailed study 
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Fig. 1. Interstellar C IV component addition: the black 
solid line shows the C iv AA1548,1551 profile of a 
CMFGEN model, and the red long dashed line is the re- 
sulting profile after including the IS absorption. A convo- 
lution has been performed to take into account the instru- 
mental resolution (1.2 A). The parameters used to model 
the interstellar component are: log(V(C IV) ga \) — 14.4, 
log(V(C IV)smc) = 14.5 and u S MC = 140 km s _1 . 

of SMC-N81. Then the value of log(iV(CIV) = 14.5 de- 
rived on average for the SMC regions in Table is chosen 
to be typical of the CIV column density in the direction 
of N81. 

With these column densities, we have created syn- 
thetic profiles of the interstellar C IV absorption lines. 
The method used was the same as that employed to esti- 
mate the column densities from the spectra of the 8 SMC 
stars mentioned above. We have then added these inter- 
stellar contributions to the CMFGEN profiles. Figure Q 
shows an example of such a correction. The typical uncer- 
tainties in the column densities translate to modifications 
of the depth of the corrected absorption profile of the order 
0.05. 

5. Detailed analysis of individual stars 

In this section, we present the results of the quantitative 
analysis of the UV spectra of our target stars. Constraints 
on the effective temperature, the luminosity, the terminal 
velocity of the wind and the mass loss rate are the main 
outputs. Secondary constraints on the slope of the veloc- 
ity field or the amount of clumping are also given. The 
method used is explained in detail for the case of star 2, 
while for the other stars the results are summarised in 
Sect. O and in Table H 
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5.1. Effective temperature 

The most reliable T c g diagnostics for O stars remain the 
photospheric Helium lines in the optical. Unfortunately, 
optical spectra of the SMC-N81 stars are not available so 
that we had to rely on the UV. Two types of indicators 
were used in this spectral range: the shape of the SED and 
the strength of several lines. 

5.1.1. UV colour index 

As O stars emit most of their luminosity in the UV, the 
shape of the spectral energy distribution at these wave- 
length is sensitive to the effective temperature just as the 
optical spectrum in the case of cooler stars. A colour in- 
dex in the UV can then allow to determine T e ff. This has 
been done on the observational side by Fanelli et al. ([1992 ) 
who have computed various spectral indices based on IUE 
spectra for different groups of stars. 

We have used the recent grid of O dwarf models by 
Martins et al. I|20f)2bjl recomputed for an SMC metallic- 
ity (Z = 1/8Zq) completed by various models at this low 
Z to derive a relation between effective temperature and a 
synthetic colour index defined by p™* 6 where F1285 (-Fisss) 
is the mean flux in an artificial 20 A wide box-shaped fil- 
ter centered on 1285 A (1585). The choice of these wave- 
lengths was a compromise between having fluxes in far 
enough wavelength ranges (to get a ratio significantly dif- 
ferent from 1) and avoiding metallicity effects (see below). 
In the UV part of the spectrum of interest to us, an in- 
crease of T c ff translates to a decrease of (whereas 
the slope of the spectrum, -F1285 — -FissSj increases). This 
is illustrated in Fig. [21 where we see that an increase of 
T GfT from 33343 K to 48529 K induces a decrease of #221 
from ~ 2.5 to ~ 1.8 while the slope increases (see the bold 
lines). This is confirmed observationally by Fanelli et al. 
fTM2l . 

Fig. |3| shows the correlation between our UV colour 
index and T c g . The determination of the UV colour index 
from the dereddened flux distribution of the SMC-N81 
stars together with this theoretical relation allows us to 
estimate T e g of the observed stars. In the case of star 2, 
we found a value of ~ 40000^. 

This method suffers from various problems. The most 
important is probably the extinction which modifies the 
slope of the SED and renders the UV colour index uncer- 
tain. The dashed line in Fig. [3]shows the position of star 
2 if its SMC extinction was increased by 0.02 (compared 
to the extinction derived from photometry): in that case, 
the estimated T c g would be ~ 34000 K. This may seem 
surprising, as naively a higher T e g could be expected for a 
higher extinction. However, this behaviour is simply due 
to the fact for a "normal" extinction law (e.g. Prevot et 
al. H984[) an increase of E(B — V) translates to an increase 
of -Fdorcd/^dcrod which corresponds to a lower T e g as dis- 
cussed above. Fig. [5] shows that a typical error of 0.02 on 
E(B-V) (i.e. an error on the flux determination) translates 




1300 1300 1400 1500 1600 1700 



MA) 

Fig. 2. UV flux distribution of two O dwarf models at 
33343 K and 48529 K. Whereas the slope Fi 285 - Fi 585 
increases with T c a the ratio ~ 12W decreases. The bold lines 
are to guide the eye and to show the variations of the slope. 

to an uncertainty of the order 3000/4000 K in T e g based 
on the UV colour index method. 

The UV SED is also shaped by many metallic (mostly 
iron) lines and thus metallicity can affect the determina- 
tion of UV magnitudes. This is clearly demonstrated in 
Fig. |21 where the effect of increasing the metallicity from 
1/8 Zq to Z strongly steepens the slope of the relation 
UV colour index - T e g. Both the slope of the continuum 
and the "line forests" are responsible for such a behaviour. 
Our choice of filters centered at 1285 A and 1585 A (where 
the density of metallic lines is weaker than in other wave- 
length ranges) was made to try to minimise the latter 
effect. 

Finally, the slope of the SED also depends on the grav- 
ity (e.g. Abbott & Hummer 1985). However, for dwarf 
stars, this dependence is weak: a test case run for a model 
with log g = 4.1 and 4.4 has shown a change of the ratio 
#2ffi of less than 2%. 

-T1585 

Given the above uncertainties, the UV colour index 
method can only give an indication of the effective tem- 
perature which needs to be confirmed by other indicators. 

5.1.2. spectral lines 

Several spectral features can be used as T e g- indicators in 
the UV. 

• O iv AA1339,1343 /O v A1371 : 

O iv AA1339,1343 is present in the spectra of most O 
dwarfs while O v A1371 appears only in stars earlier than 
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Table 3. Effective temperature estimates. 



O 



A 


A 


i.ii 




- 




A 


- 


- 

▲ 






star 3 








star 2 E(B-V) auc + 0.02 
star 1 1 






A A 


A star 1 






A 


star? 






A A 


▲ 


I i_ 




1,11, 


▲ 

A ▲ 

i , i I , , i 



3x10* 3.5x10* 4x10* 4.5xl0 4 5x10* 

T et , (K) 

Fig. 3. Effective temperature indicator: colour index 
method. The colour index is define by ~^ 1285 where F1285 
(-F1585) is the mean flux at 1285 A (1585 A). Filled (open) 
triangles are for models with Z = 1/8Z@ (Z = Zq). 
Horizontal lines indicate the colour index for N81 stars 
which, by comparisons with the theoretical values, give 
an estimate of the effective temperature. Note that this 
method is strongly metallicity-dependent. The dashed line 
show the position of star 2 for a SMC extinction increased 
by 0.02. 



06 (Walborn et al. 11995(1 so that their presence/ absence 
and relative strength is a good T e g estimator. Fig. 01 
shows comparisons between models of different T c g and 
the observed spectrum of star 2. An increase of T e g in- 
creases the strength of both lines which are reproduced 
for T cS > 40000 K. 

O v A1371 alone has been used by de Koter et al. 
(1998) as a T c g indicator for early O stars. Nonetheless, 
as mentioned by these authors, this line is always pre- 
dicted too strong compared to observations above ~ 40000 
K. Below this limit, the situation seems better. Recently, 
Bouret et al. PQ03jl have claimed that clumping could help 
to solve this well known problem (see Schaerer et al. 119961 
Pauldrach et al. I2001[) : the line is weaker in a clumped 
model than in an homogeneous wind. As O v A1371 (to- 
gether with O iv AA1339,1343 ) depends also on M, our 
T e ff estimate relies on the relative strength of both lines 
in homogeneous winds. This estimate must be confirmed 
by stronger indicators. 

• Fe IV/Fe V: 

Iron lines forests exist all over the UV range. In partic- 
ular, Fe IV lines are present between 1600 A and 1640 
A, while Fe V lines are found between 1430 A and 1480 



Star 


Estimator 


T off (K) 


1 


UV colour index 


37000 




OIV/OV 


36000 




Fe IV/Fe V 


> 38500 




N V 






C III 


38500 


2 


UV colour index 


37500 




O IV/O V 


> 40000 




Fe IV/Fe V 


40000 




N V 


> 35000 




C III 


42000 


3 


UV colour index 


34000 




O IV/O V 


< 37500 




Fe IV/Fe V 


36000 




N V 


< 37500 




C III 


36000 


11 


UV colour index 


35000 




O IV/O V 


37000 




Fe IV/Fe V 


> 37000 




N V 


< 37000 




C III 


37000 



Adopted T cff (K) 



38500 



40000 



36000 



37000 



A. The iron ionisation increases with T g so that Fe IV 
lines weaken relatively to Fe V lines between 35000 and 
42000 K. As Fe V is the dominant ionisation state of Fe 
in this temperature range, Fe V lines are saturated and 
little affected by an increase of T e g whereas Fe IV lines 
weakens. This is shown in Fig. 0] For star 2, a value of 
at least 40000 K is necessary to reproduce the observed 
spectrum. 

The determination of T e g from this line ratio can be 
hampered mainly by two effects: 

- the iron abundance (and more generally metallicity) 
can change the strength of the Fe absorption. This effect is 
twofold. First, increasing the iron abundance will immedi- 
ately increase the absorption of all Fe ions, although differ- 
ently depending on the position of the lines on the curve of 
growth. Consequently, the ratio of lines from two succes- 
sive ionisation states will be modified. Second, increasing 
the iron abundance will strengthen the line-blanketing ef- 
fect and thus will increase the local temperature in the line 
formation region. The ionisation will be increased, lead- 
ing to a higher T e g estimate. The effects of metallicity due 
to line-blanketing on the effective temperature of O stars 
have been estimated by Martins et al. (2002b|) and turn 
out to be of the order 1000 to 2000 K depending on the 
spectral type for metallicities ranging from solar to 1/8 
solar. We have run test models for a global metallicity of 
Z = 1/5Zq. Its main effect is to strengthen the Fe IV 
lines, leaving the Fe V lines unchanged as they are al- 
most saturated in the temperature range of interest here. 
Quantitatively, this change of Z is equivalent to a decrease 
of T e ff by ~ 2000 K as regards the iron lines behaviour. 
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- the so-called turbulent velocity "artificially" 
strengthens the absorption profiles as shown in Fig [121 
where «turb increases from 5 to 15 km s _1 . Fe IV and 
Fe V lines deepen differently when i>t U rb increases, so 
that their relative strength remains turbulent dependent. 
Nonetheless, the effect is smaller than the effect on indi- 
vidual lines so that a reasonable estimate of T e g can be 
drawn from the study of the relative strength of these iron 
forests. Note that if the effective temperature is known, 
the iron lines can be used to determine the turbulent ve- 
locity (see Sect. 15. 511 . 

• C in AA1426, 1428 : 



As noted by Walborn et al. (j!995f) in their IUE atlas of O 
star spectra, this blend of C III lines strengthens towards 
later types. Fig.0]shows that in the case of star 2 a value of 
at least 42000 K is required to fit the observed spectrum. 
Changing slightly the carbon content can affect this deter- 
mination. Quantitatively, a reduction of the C abundance 
from 1/8 C Q to 1/10 C© (as discussed in section l5".2.2f) is 
equivalent to an increase of T c g by ~ 1500 K. 

• N v A1238,1242 : 

This resonance line is known to be the strongest around 
spectral type 04 (Smith-Neubig & Bruhweiler 1997) and 
to weaken at later spectral types or equivalently when T c g 
decreases. Quantitatively, a T e g > 35000 K is required to 
account for the N v A1238,1242 absorption profile in star 
2. Nonetheless, the strong mass loss dependence of this 
line makes it a poor and only secondary T e g indicator: a 
high mass loss associated with a relatively low T e g can 
mimic the NV profile of a star with a lower mass loss rate 
but a higher T c g. 

Table summarises the results of these T e g estimates 
which all point to an effective temperature of the order 
40000 K for star 2. A real dispersion exists and is mostly 
due to the multi-parameter dependence of most of the 
indicators used. 



5.1.3. Estimation of uncertainty on T e g 

Table shows that for a given star, the dispersion in the 
T e g estimators is of the order ± 2500 K in the worst cases. 
A possible additional source of uncertainty comes from 
the rectification of the spectra. We estimate it to be lower 
than 0.05 in term of normalised flux, which from Fig. 0] 
(especially the FeV plot) can lead to an error on T e g of 
the order ± 1500 K. 

The uncertainties due to variations of the global metal- 
licity and carbon content have been estimated above and 
are of the order 2000 K. 

To estimate the efficiency of our T e g determination 
method, we have applied it to a star for which optical 
and UV analysis have given strong constraints on T e g. 
For this purpose, lOLac was chosen (see section l7~2f . The 
UV colour index method indicate a value of 39000 K, while 
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Fig. 4. Spectroscopic T e g indicators. The black solid line 
is the observed spectra of star 2. Coloured lines are four 
different models with T c g 42000 K (red, dot-long dashed 
line), 40000 K (magenta, short dashed line), 37500 K 
(green, long dashed line), 35000 K (blue, dotted line). The 
mass loss rate is fixed at a constant value of 10 -85 M Q 
yr _1 . See text for discussion. 



the spectral lines point to values of the order 35000/36000 
K. As the accepted T e g for lOLac is close to 36000 K, the 
error on the T c g estimate is not more than ± 3000 K. This 
example also shows that the spectral lines method is more 
accurate than the colour index method. For that reason 
we have given more weight to the T c g estimates based on 
the line method in our final estimate (as can be shown in 
Table El). 

From the above discussion, we conclude that on 
average, the uncertainty on our T e g determination is of 
the order ± 3000 K. 



5.2. Mass loss rate 

The determination of the mass loss rate remains one of the 
main goals of this study. As already noted, the spectra of 
N81 stars show very weak wind features with no emission. 
As emission, contrary to absorption, is entirely formed 
in the wind, our constraints on M have been set by the 
requirement that no emission is produced in the model 
spectra (as observed in the STIS spectra). Thus, we have 
derived only upper limits on the mass loss rates as models 
with M below this limit never produce emission. The main 
mass loss indicators in the UV are the N v A1238,1242 , 
C iv AA1548,1551 and O V A1371 lines. 
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5.2.1. Primary determinations: 

• O v A1371 : 

O v A1371 develops a P-Cygni profile in the earliest super- 
giants while only an absorption is seen in dwarfs (Walborn 
et al. 119950 . Fig. |S]reveals that the profile deepens when 
M > 10 -8 M Q yr -1 but remains roughly unchanged be- 
low this value. Fitting the observed spectrum of star 2 
requires a mass loss rate lower than 10 -8 M Q yr -1 . As 
mentioned in the previous section, O v A1371 is sensitive 
to T e g so that a degeneracy M /T e g exists. X-rays can also 
affect this line by increasing the Oxygen ionisation to pro- 
duce O VI (observed in the EUV) by Auger ionisation of 
O IV. 

• N v A1238,1242 : 

This line is present in dwarfs of spectral type earlier than 
08 (Walborn et al. ITMSI Smith-Neubig & Bruhweiler 
I1997fl and shows a strong P-Cygni profile in early dwarfs. 
It is mostly formed in the wind so that it is one of the 
best mass loss rate indicators of the UV part of the spec- 
trum. This is illustrated in Fig. Elwhere we see the profile 
changing from a well developed P-Cygni shape to a weak 
absorption when M is reduced by two orders of magni- 
tudes. A mass loss rate lower than 10 -9 M Q yr -1 is nec- 
essary to produce no emission. However in that case, the 
absorption profile is too weak. Increasing M by a factor of 
3 improves the fit of this absorption part but induces an 
emission (see Sect. for a discussion of the shape of these 
wind profiles). Nonetheless, as emission is only produced 
in the wind whereas absorption can originate both in the 
wind and in the photosphere, we prefer adopt M = 10 -9 
Mq yr -1 as a reasonable upper limit. 

N v A1238,1242 is also T e g sensitive as seen in Sect. 
15. 1.21 so that any error on T e g can lead to an error on M. 
Another problem comes from the X-rays (supposed to be 
created by shocks in the outer wind, e.g. Owocki et al. 
1988, Feldmeier et al. H997fl which can increase the ionisa- 
tion of Nitrogen, thus leading to modifications of the N v 
A1238,1242 profile. A very accurate determination of M 
based on N v A1238,1242 would then require the inclusion 
of X-rays in the models. However, tests have revealed that 
the inclusion of X-rays does not modify our conclusions. 
It is also unlikely that X-ray emission is important since 
any emission in the wind of the SMC-N81 stars seems to 
be very weak as shown by the UV spectra. 

• C iv AA1548,1551 : 
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Fig. 5. Mass loss rate indicators: the black solid line is the 
observed profile while coloured lines are taken from models 
with T eff =40000 K and various M (-log(M)=7 (red, short 
dashed line), 7.5 (yellow, long dashed line), 8 (magenta, 
dot-short dashed line), 8.5 (green, dot-long dashed line), 
9 (blue, short dashed - long dashed line)). See text for 
discussion. 



the emission disappears. But this absorption is mainly 
photospheric and the wind part turns out to be too 
weak compared to the observed spectra. A change of T e g 
does not lead to any improvement (a lower T e g implies 
an emission in C iv AA1548,1551 and a too weak N v 
A1238,1242 line, and a higher T e g weakens further the 
C IV AA1548,1551 line). 



5.2.2. Improving the M determinations: effects of (3, 
clumping, adiabatic cooling and abundances 

A difficulty to fit the UV lines (especially C iv 
AA1548,1551 and N v A1238,1242 ) is to produce a sig- 
nificant absorption extending up to velocities > 1800 km 
s _1 without any emission. In the following, we concentrate 
on C IV AA1548,1551 as the discrepancy is highest for this 
line, but the discussion applies equally to N v A1238,1242 



This line is the other strong UV mass loss indicator 
of O dwarfs. It is seen in all O dwarfs with a strength 
increasing towards early types (Walborn et al. 11995(1 . 
In Fig. [5] we see that it shifts from a P-Cygni profile 
to an absorption profile when M decreases. The absence 
of emission, required by the observation, is obtained for 
M < IQ- 7 - 5 M yr- 1 . Contrary to N v A1238,1242 , 
a significant absorption profile can remain even when 



The absence or weakness of absorption at high velocity 
in the models showing no emission results from a lack of 
absorbers (i.e. C IV ions) which may come either from a 
too low density, and then from a too low mass loss rate, or 
from a too high ionisation of carbon (CV being the domi- 
nant ionisation state in the wind) . The first possibility can 
be ruled out because a higher mass loss rate will produce a 
non observed emission. Relying on the second hypothesis, 
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we have sought for mechanisms that could lead to a reduc- 
tion of the ionisation in the wind. This can be achieved 
if recombination rates are increased, i.e. if the density is 
higher (recombination scales as p 2 ). As the density is given 
by p = 47T ^2 v f where f is the filling factor (/ =1 in an ho- 
mogeneous atmosphere) , a higher density at a given radius 
can be obtained by either a stronger mass loss rate (which 
is excluded), or a clumped wind (/ < 1) or a lower ve- 
locity which, as Voo is fixed, implies a softer slope of the 
velocity field (the (3 parameter). Recombinations can also 
be increased when the temperature in the outer wind is 
reduced. Adiabatic cooling may induce such a reduction 
as in these low density winds it can become an important 
cooling process. Last, abundances can of course modify the 
strength of the wind profiles. In the following, we present 
investigations of the influence of these various parameters 
on the line profiles and the mass loss rates determinations. 

• (3 effects: 
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P is usually determined through the shape of Hydrogen 
emission lines in the optical range (e.g. Hillier et al. 2003). 
However, as optical spectra are not available and as the 
optical lines have probably absorption profiles due to the 
weakness of the winds, we have no constraints on this pa- 
rameter. We have then run test models with (3 — 2.0 (the 
default value in all our computations being 0.8). The re- 
sults are displayed in Fig. C iv AA1548,1551 is not mod- 
ified since it is almost purely photosphcric in the models 
shown here, while N v A1238,1242 shows narrower and 
stronger absorption and emission when /? is higher. This 
behaviour can be explained in term of size of the inter- 
action region which is the region in which, due to the 
wind velocity induced doppler shift, a photon can interact 
with a given line. The size of this region is proportional to 
the Sobolev length (see Lamers & Cassinelli 1999) which 
scales as (dv/dr) -1 in the inner wind, so that a higher /3 
will lead to larger Sobolev length (the acceleration being 
smaller) . The interaction region is then wider which leads 
to a stronger absorption or emission in the center of the 
line. In the outer atmosphere, the radial Sobolev length is 
proportional to f3~ 1 r 2 /v (for a (3 velocity law) and smaller 
for higher (3 so that the emission/absorption are reduced 
at high velocity 1 . Illustrations of the dependence of the 
size of the interaction region can be seen in Fig. 8.6 of 
Lamers & Cassinelli l|1999[) . 

In our case, low values of (3 seem to be preferred as the 
observed N v A1238,1242 profile does not show a double 
absorption line. However, (3 — 2.0 reproduces better the 
blue edge of the absorption profile. As f3 — 0.8 is closer 
to the predictions of the radiation driven wind theory, we 
adopt this value as typical of the SMC-N81 stars. But 



1 When (3 increases, the transverse Sobolev length, which is 
proportional to (v/r)~ l , increases too. However, the decrease 
of the radial Sobolev length is higher than the increase of the 
transverse Sobolev length so that globally, the size of the in- 
teraction region is reduced. 



Fig. 6. Influence of the slope of the velocity field ((3) on 
the mass loss rate diagnostic lines. Models with (3 = 2.0 
(blue short dashed line) and (3 — 0.8 (red long dashed line) 
are compared to the observed profiles of star 2 (black solid 
line). Models are for T eS = 40000 K, M = 10~ 8 - 5 M© yr" 1 
and V sim = 300 km s -1 . The C iv AA1548,1551 line is 
almost unchanged while the N v A1238,1242 line shows 
narrower absorption and emission when j3 is higher. 

whatever the exact value of /?, the mass loss rate deter- 
mination is not strongly modified as the level of emission 
remains roughly the same (see Fig.^J. 

• clumping effects: 

As already mentioned, there are indications that the wind 
of O stars are clumped (Crowther et al. 120021 Hillier et 
al. IM)3l Bouret et al. IMBl Repolust et al. VMty . al- 
though there exist no quantitative constraints. The ef- 
fect of inhomogencous winds are twofold: first, due to 
the presence of over densities, the emission of density- 
sensitive lines is strengthened; second, the higher density 
in clumps increases the recombination so that the ionisa- 
tion is reduced. The competition between the two effects 
can lead to either stronger or weaker emission. To inves- 
tigate deeper the effect of clumping on the wind profiles, 
we have run test models with = 0.01 (recall that the 

filling factor f is given by / = + (1 - /oo)e^). The 
main effects on the N v A1238.1242 and C iv AA1548,1551 
lines are shown in Fig. [7| We see that the strength of N v 
A1238,1242 is reduced and that the C iv AA1548,1551 ab- 
sorption is increased in the outer part of the atmosphere. 
This is attributed to a reduced ionisation in the outer wind 
(N V recombines to N IV and C V to C IV), which im- 
proves the fits. Hence, there is no doubt that the inclusion 
of clumping is crucial to reproduce the observed features. 



10 



Fabrice Martins et al.: Young massive stars in SMC N81 






Fig. 7. Influence of clumping on the mass loss rate diag- 
nostic lines. Models with a volume filling factor of 0.01 
(blue long dashed line) and 1.0 (red short dashed line, no 
clumping) are compared to the observed profiles of star 
2 (black solid line). Models are for T cff = 40000 K and 
M = 10~ 8 - 5 M yr" 1 and V ami = 300 km s _1 . The 
strength of N v A1238,1242 is reduced and the absorp- 
tion of C iv AA1548,1551 is increased when clumping is 
included. 



Fig. 8. Effect of the inclusion of adiabatic in model com- 
putations on the temperature structure (top) and the C IV 
AA1548,1551 profile (bottom; blue long dashed line: model 
with adiabatic cooling; red short dashed line: model with- 
out adiabatic cooling). While the temperature is strongly 
reduced in the outer parts, the C IV AA1548,1551 profile 
is only slightly affected. Other wind lines do not show any 
change. 



Whether this is a proof of the inhomogeneity of O star 
winds or just a trick to simulate the correct ionisation is 
not clear, but this parameter trans out to be an impor- 
tant ingredient of the modeling. As regards the mass loss 
determination, Fig. shows that the emission is slightly 
reduced when clumping is included so that the upper limit 
on M must be slightly increased. 

• adiabatic cooling: 

In low density winds, due to the reduction of any cooling 
processes based on atomic mechanisms, adiabatic colling 
is expected to be important to set the temperature struc- 
ture (e.g. Drew ir^BSfl . Fig.[S]demonstrates that it is indeed 
the case: a model with adiabatic cooling (long dashed line) 
shows a strong drop of the temperature in the outer wind. 
However, as in this part of the atmosphere populations 
are mostly governed by radiative processes, the influence 
on the line profiles is reduced: the C iv AA1548,1551 line 
shows a slightly enhanced absorption which improves the 
fit but remains marginal. We conclude that adiabatic cool- 
ing is not a crucial parameter as regards the fit of UV wind 
lines. 

• abundances: 

Metallicity in our models has been chosen to be 1/8 solar, 
and the individual abundances have simply been scaled 



according to this global metallicity. A better assumption 
would have been to take abundances typical of the SMC 
molecular clouds since the N81 stars are young and 
their atmospheres probably not contaminated by stellar 
nucleosynthesis products. Such initial abundances can be 
obtained by studies of either main sequence stars which 
have not yet experience internal mixing, or HII regions. 
However, few of the former have been undertaken and 
are often uncertain (Venn 1999 and references therein), 
and studies of nebular abundances can lead to underes- 
timates due to the depletion of material on dust grains. 
Nonetheless, reasonable agreement between these two 
types of determinations can sometimes be obtained (see 
Venn 1999 for a discussion) and point to the following 
values: C/C Q = 1 /10, N /N Q = 1/20 a nd O/Oq = 1/5 
(Venn IT39T41 Heap 1271(131 Vermeij EES)- Adopting these 
abundances in our computations leads to weaker N v 
A1238,1242 and C iv AA1548,1551 profiles for a given 
M. This means that our upper limit on the mass loss 
rate necessary to remove emission in the wind lines has 
to be increased. Recently, Asplund 2003] has revised the 
solar N abundance. With his new value, the SMC-N81 N 
abundance derived by Vermeij et al. (2002) is 1/30 solar, 
meaning that the upper limits on M has to be increased 
again a little more. 
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From the above discussion it results that if the CNO 
abundances are reduced and set to more realistic values 
and if clumping is included, the upper limit on M has to 
be increased by ~ 1.0 dex (the effect of abundances being 
dominant). Clumping also helps to get better shapes of 
the wind lines. (3 and adiabatic cooling have almost no 
influence on the M determination. As a consequence, we 
adopt an upper limit on the mass loss rate of star 2 of 
1CT 8 - M yr- 1 . 

5.3. Reliability of the mass loss rate determination 

In view of the low values of M derived, we may won- 
der whether our determination is not hampered by any 
modeling problem. In particular, what is really determined 
through fits of UV wind lines is the product of the mass 
loss rate time the ionisation fraction (qi) of the absorb- 
ing/emitting ion. Any problem with the prediction of these 
ionisation fractions would translate in an error on M. To 
investigate this point, two kinds of tests have been pur- 
sued. 

• H a vs UV mass loss determination: 

H a is much less sensitive to the model predictions con- 
cerning the ionisation than UV resonance lines, so that we 
have compared the values of M derived from the fit of UV 
lines on the one hand and H a on the other hand. For 
this aim we have chosen as a test case the star HD2 17086 
(07Vn) for which constraints on the mass loss rate (Puis 
at al. H99ril Repolust et al. ??) and the ionisation fractions 
of several ions (Lamers et al. 1199 91 exist. The low density 
of the wind (see Lamers et al. I1999H and the effective tem- 
perature of the order 37000 K make this star similar to 
the SMC-N81 stars, although the wind is less weak. We 
have computed various models to fit simultaneously the 
H a profile and the UV resonance lines. 

Fig.[5]shows the results of the fits of the wind sensitive 
lines. Two types of conclusions can be drawn depending 
on the value of the (3 parameter (slope of the wind velocity 
field): 

- if we adopt (3 — 0.8 as derived by Repolust et al. 
the H a profile is best fitted for M = 10~ 6 - 4 M 



yr (dot dashed line) but in that case the UV lines (espe- 
cially N iv A1718 ) are too strong. A value of M of 10~ 7 ' 2 
M© yr -1 (long dashed line) has to be adopted to improve 
the fit of these UV lines, but now the H a absorption is too 
strong. Mass loss rates derived from the UV seem then to 
be lower by a factor 6 compared to the H a determina- 
tion. Note that the predictions of Vink et al. I)2001|l gives 
M = 10~ 61 Mq yr" 1 for HD217086. 

- when we increase (3 up to 1.7, a mass loss rate of 
10~ 70 M yr -1 gives a reasonable - although not perfect 
- agreement between observations and models for both 
H a and the UV lines (dotted line on Fig. |SJ. 

While the conclusions of the case (3 — 0.8 point to 
a problem with the ionisation fractions predicted by the 
CMFGEN models, the case /3 = 1.7 indicate that this 



problem partly disappears when (3 is increased. This does 
not necessarily mean that (3 — 1.7 (the determination of 
this parameter in low density winds is difficult, see e.g. 
Puis et al. H996f> . but it is a mean to produce an ionisation 
throughout the wind which leads to an overall good fit of 
the observed profile. 2 Practically, as we do not have any 
constraint on (3 (see Sect. 15.2.2(1 we can conclude that in 
the worst case M derived from UV lines is underestimated 
by a factor 6, at least if we take the mass loss rate from 
H a as the correct one. In fact, if we take the case f3 = 1.7 
as representative of the real conditions in the wind, the 
error we make when we only look at the UV spectrum (as 
for our SMC-N81 stars) and adopt (3 = 0.8 is of 0.2 dex, 
or less than a factor 2. 3 Moreover, in their recent study of 
SMC dwarfs with weak winds, Bouret et al. H200ri[) have 
been able to fit simultaneously lines of different ionisation 
stages of the same element with CMFGEN models which 
indicates that the wind ionisation was correctly predicted. 

• SEI method: 



We have applied the SEI method (Lamers et al ll987|l to 
our SMC N81 star 2. This method leads to the determi- 
nation of M x (jj. Basically, the SEI method solves the ra- 
diative transfer in an expanding atmosphere for which the 
source functions are calculated in the Sobolev approxima- 
tion. The main input parameters are the velocity field and 
a function giving the optical depth of the line as a function 
of velocity. The main output is the value of M x q\. The 
best fit to the C iv AA1548,1551 line of star 2 is shown in 
Fig. ED The corresponding value of log(M x q\) is -9.68. 
We also show in this figure the influence of the underly- 
ing photospheric absorption: as expected, the low-velocity 
part of the C IV AA1548,1551 absorption profile is sensi- 
tive to this photospheric absorption. As the interstellar 
contamination renders the true stellar absorption uncer- 
tain (see above), we do not give too much weight to this 
part of the profile which can always be fitted by tuning 
the photospheric absorption. The important point is that 
the wind profile (0.4 < Av/voo < 1.0) is insensitive to the 
amount of photospheric absorption and only depends on 
the wind parameters. Our estimation of M x qi with the 
SEI method can then be regarded as reliable. 

In order to compare this value to that deduced from 
our model giving the best fit for star 2, we have com- 
puted the ionisation fraction of C IV using the definition 
of Lamers et al. (1999): 



5i 



Jxq nj{x)dx 
C n E (x)dx 



(1) 



2 As mentioned in Sect. I5,231 clunrping can also modify the 
ionisation in the wind. However, several tests have shown that 
it was not possible to fit simultaneously H a and UV lines with 
clumped models. 

3 Compared to the value given by Repolust et al. f 2003 1 - 
M = 10" 6 64 M yr" 1 -, our best determination (M = 10" 7 
Mq yr~ , (3 = 1.7) is a factor ~ 2 lower 
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Fig. 9. H a versus UV lines M determinations. The solid 
line is the observed profiles of wind sensitive lines of 
HD2 17086. Other lines are CMFGEN models at 37000 
K with /3 = 1.0 and M = 10~ 6 - 2 M Q yr" 1 (dot-dashed 
line), P = 1.0 and M = 10~ 7 - 2 M Q yr" 1 (long dashed 
line) and (3 = 2.0 and M = 10~ 6 - 8 M Q yr" 1 (dotted line). 
A rotational velocity of 290 km s _1 has been adopted. For 
j3 = 1.0 values of M different by one order of magnitudes 
are required to fit H a on the one hand and the UV lines 
on the other hand, while for (3 = 2.0 a reasonable fit of all 
lines is achieved. See text for discussion. 



where nj is the population of the absorbing ion, tie the 
population of the element, xq (xi) the lower (higher) inte- 
gration limit (x being r/R+). We found logqc iv = —1-36 
which, together with a mass loss rate of 10 -85 M yr _1 
gives log(Af<7i) = -9.86, in reasonable agreement with the 
SEI result. This shows that the value of M x q\ obtained 
by fitting the CMFGEN model profiles to the observations 
is correct. Thus it follows that if M is underestimated, the 
C IV ionisation fraction is overestimated. To investigate 
this point, we have compared our ionisation fractions to 
those derived by Lamers et al. I|1999[) . For stars of different 
stellar and wind properties, they found a mean qc rv of 
the order 10 -25 which is ~ a factor of 10 lower than the 
CMFGEN value. An error by such a factor on g; translates 
to an underestimate of M by the same factor. This con- 
firms the result of the previous section where a possible 
underestimate of M by a factor 6 was highlighted. 

An interesting comment to make is that if we simply 
assume that the derived mass loss rate is correct, the 
mean density in the wind (as defined by Lamers et al. 
1999) is of the order io -16 — -17 g cm~ 3 which is well 
below the lowest density probed by the Lamers et al. 
sample (see their Fig. 3). As their analysis indicates an 



increase of the C IV ionisation fraction with decreasing 
mean density, it is conceivable that for very weak winds, 
qc iv m ay reach values of the order 0.1, which would 
reconcile the CMFGEN and Lamers et al. ionisation 
fractions. Moreover, it should be noted that the study of 
Lamers et al. includes radio and H a mass loss rates which 
did not take clumping into account. Since the inclusion 
of clumping leads to lower M (e.g. Hillier et al. 2005|), it 
is conceivable that some of the Lamers et al. ionisation 
fractions are underestimated. 

The main conclusion of these two types of studies 
(H Q /UV and SEI) is that there are indications that the 
ionisation fractions predicted by CMFGEN may be wrong, 
but by no more then a factor of ~ 10 (leading in that case 
to an underestimation of M by the same factor). However, 
in both studies we have also found means to explain the 
observations with the CMFGEN predictions (changing the 
value of (3 in the study of HD217086, extrapolating the 
trend qc iv - mean density in the SEI study) . It is thus not 
clear if the ionisation fraction predictions of the CMFGEN 
models are erroneous or not, especially given that we are in 
a range of parameters never explored before (weak winds) . 
However, if these predictions are wrong, they do not qual- 
itatively modify the conclusion concerning the weakness 
of the winds (see Sect. 0). 

5.4. Terminal velocity 

The terminal velocity of the winds of O stars is usually de- 
rived from the blueward extension of UV resonance lines. 
Here, due to the weakness of the outer wind density, the 
absorption may not extend up to so that with the 
above method one can only derive lower limits for the ter- 
minal velocities. As shown in Table 0] these limits can 
even be lower than the escape velocity (which is of the 
order 1100 km s _1 ), reinforcing the fact that the absorp- 
tion probably does not extend up to v^. Moreover, the 
relatively low signal to noise ratio of our spectra coupled 
to the uncertainty in the flux normalisation introduces 
an uncertainty in the exact position of the most blue- 
shifted absorption. Figure ITD shows the N v A1238,1242 
and C iv AA1548,1551 profiles of models with Voo = 1500 
and 1800 km s _1 compared to the observed line. The N v 
A1238,1242 profile is not affected by the change of Voo 
while the C iv AA1548,1551 profile is hardly modified. In 
view of this result and as the fit with the SEI method re- 
quires a value of 1800 km s _1 for the terminal velocity, we 
adopt this value as representative for star 2. The values 
for the other stars are given in Table 0] 

5.5. Turbulent velocity / Rotational velocity 

The turbulent velocity in O stars is thought to increase 
from values of a few km s _1 near the photosphere to 
~ 10% of the terminal velocity in the outer atmosphere. 
As mentioned in Sect. 15. ll iron lines are sensitive to i>t U rb 



Fabrice Martins et al.: Young massive stars in SMC N81 



13 



1.2 
1 




i 1 i i 


1 


■/.■ 


.,11 

■ 






v U 




u 




0.8 








if 

■i 

ii 




0.6 






\\ 

\ ■ 






0.4 






; V" 


! J 




0.2 




1 




■ 


i 



-1 1 

<5(v)/v 

Fig. 10. Fit of C iv AA1548,1551 line with SEI method. 
The solid line is the observed profile and the dashed line 
is the best fit. The parameters for the velocity law are 
f3 = 1.0 and «oo = 1800 km s _1 . A photospheric compo- 
nent has been used. The optical depth has the following 
velocity dependence: t(w) = T tot (l — (w) 1 ^)(w) - 2 with 
w = v/voo and T tot = J^ 01 t(w). The derived value of 
M x q; is 10 -9 68 . Dotted lines show the influence of an 
increased or reduced photospheric component: the wind 
part of the absorption is not affected by this component. 
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so that once T c g is known, they can be used to determine 
the turbulent velocity. Fig. 1121 shows the comparison of 
iron lines of a model with T c g = 40000 K but different 
fturb- For each model, spectra convolved with rotational 
velocity of 200, 250 and 300 km s" 1 are shown. The best fit 
is obtained for turbulent velocity of 5 km s _1 . For higher 
values, the Fe V lines deepens too much compared to the 
observed spectrum. The Fe IV lines behave similarly but 
also weakly. Tests have been run with a turbulent velocity 
varying from 5 km s _1 near the photosphere up to 100 km 
s _1 in the outer wind and have revealed very little changes 
in the blue part of wind profiles. This is explained by the 
weakness of the wind in which absorption is likely to take 
place up to velocities lower than the terminal velocity of 
the wind. Fig. ^] also reveals that the rotational veloc- 
ity of the star is of the order 300 km s _1 . Indeed, lower 
velocities lead to narrow and deep profiles not observed. 
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Fig. 12. Effect of microturbulence and rotation on the 
iron spectrum. The three left panels show Fe V lines while 
the Fe IV lines are shown on the right panels. The tur- 
bulent velocity decreases from 15 to 5 km s _1 from top 
to bottom. Each model is convolved with a rotational ve- 
locity of 100 (dot - long dashed line), 200 (short dashed 
line) and 300 km s _1 (long dashed line). The lower panels 
(«turb = 5 km s _1 ) give the best fit. 
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5.6. Luminosity 

The luminosity of the stars has been estimated from the 
absolute visual magnitude My and a bolometric correction 
calculated for the estimated T e g of the star according to 

log-^- = -0A(M v + BC - M*° l ) (2) 

Lq 

where M| oi = 4.75 (Allen 1976}. M v was derived in 
paper I from the observed visual magnitude and the esti- 
mated extinction. BC, which is essentially model indepen- 
dent when calculated as a function of T c g , is derived from 
the relation of Vacca et al. I|1996[) 

BC(T eS ) = 27.66 - 6.84 x log T cS (3) 

The uncertainty is 0.01 for the observed visual magni- 
tude (paper I), 0.05 for the extinction (from an uncertainty 
of 0.01 on E(B — V)), 0.025 for the distance modulus (di 
Benedetto I1997|) and finally 0.25 for BC for a typical er- 
ror on T e s of 3000 K. On average, the uncertainty on L is 
therefore of the order 0.12 dex. 

5. 7. Summary and results for other stars 

Stellar and wind parameters for 3 other stars of SMC-N81 
have been determined with a similar analysis. Here we 
summarise the results of this analysis and give the results 
in Tabled The best fits are shown in the Appendix. For 
the remaining stars presented in paper I, no constraints 
have been derived due to the poor quality of the spectra. 
For all the stars, (3 = 0.8 and Si, S and Fe abundances 
equal to 1/8 solar have been chosen. The masses have 
been derived from the HR diagram presented in Fig. 

5.7.1. star 2 

The main results have been given in the previous sections. 
In contrast to all other stars, we stress that the fits of the 
C iv AA1548,1551 and N v A1238,1242 lines for star 2 
are improved when the abundances are taken 1/8 solar. 
Choosing the lower values indicated by the nebular anal- 
ysis (see Sect. 15.2. 2|l lead to too weak absorptions even 
when clumping is included. The best fits are achieved with 
a filling factor at the top of the atmosphere of 0.01 4 . 

5.7.2. star 1 

For this star, T e g is found to be ~ 38500 K. The upper 
limit on the mass loss rate is 10~ 85 M yr _1 and i) m is at 
least 1500 km s _1 . Clumping is necessary to improve the 
fit of the C iv AA1548,1551 line (/„ = 0.01). A reasonable 
fit is achieved with the nebular CNO abundances. 

4 Recall that in the present parameterisation of clumping in 
CMFGEN, this means that / goes from 1 at the photosphere 
to 0.01 when v = Voo. 



5.7.3. star 3 

Star 3 is the coolest of the 4 stars studied (T e g = 36000 
K). The determination of the mass loss rate is very diffi- 
cult because the N v A1238,1242 line is almost absent due 
to the low T e ff and the C iv AA1548,1551 line is almost 
entirely photospheric. Hence, secondary indicators such 
as O iv AA1339,1343 have been used estimate M (which 
is found to be lower than 10 -85 M Q yr _1 ). We want to 
stress that this determination is probably the most un- 
certain of our sample, especially since the the observed 
spectrum is noisy. The very small value for the terminal 
velocity reflects the low density of the wind: the number 
of absorbant in the wind is so low that the wind lines are 
essentially absent. The nebular CNO abundances give the 
best fits. 

5.7.4. star 11 

We have estimated an effective temperature of 37000 K for 
this star. An upper limit of 10~ 95 M© yr~* is estimated 
from the N v A1238,1242 and C iv AA1548,1551 lines. 
The terminal velocity we derive (600 km s~ x ) is again a 
lower limit. The nebular CNO abundances give reasonable 
fits and clumping is not necessary. 

6. Nebular and stellar properties 

In this section we first go back to the nebular properties 
of N81 and then investigate the evolutionary status of the 
individual stars together with the consequences of their 
weak winds. 

6.1. Nebular properties / Ionising fluxes 

A meaningful way of testing our results concerning the 
stellar properties of the N81 individual components is to 
compare them to the integrated properties of the cluster. 

First, Heydari-Malayeri et al. (1999) derived a mean 
extinction of Ay = 0.40 from the observed H a /Hp ratio 
in the nebula, and they use this value to correct the ob- 
served flux in Hp and to estimate the number of Lyman 
continuum photons - Qq - emitted by the N81 stars (under 
the assumption that the HII region is ionisation bounded). 
They find Q = 1.36 x 10 49 photons s" 1 . Fr om this and 
the calibration of Vacca et al. H1996|) . they conclude that 
a single main sequence star of spectral type 06.5 or 07 
can lead to such an ionising flux. We have estimated the 
total amount of ionising photons released by the N81 stars 
from the SEDs of the models giving the best UV fit. As the 
N81 stars studied here are the most luminous and hottest 
of the region, they are likely to provide essentially all the 
ionising flux. We found from the models Qo = 1.64 x 10 49 
photons s" 1 in good agreement with the values derived 
from the nebular properties. The lower value of the esti- 
mate of Heydari-Malayeri et al. ( 199HJ) can be partly ex- 
plained by the fact that the HII region may be density 
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Table 4. Summary of the stellar and wind properties derived for the N81 stars. The spectral types are estimated 
from the optical spectra of the models giving the best UV fits. As we have only lower limits on the terminal velocities, 
we have adopted Yoo/^esc = 2.6 to compute both the modified wind momenta and the theoretical mass loss rates of 
Vink et al. I|2001[l . The Si, S anf Fe abundances are 1/8 solar and n(He)/n(H) is 0.1. The gravity adopted for our 
computation was log g = 4.0 as it is typical of dwarf stars (Vacca et al. H996[l and as we have no diagnostics to estimate 
the value of this parameter. 
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bounded so that a part of the ionising flux may escape 
the cavity. 

Second, we have integrated optical spectra of the N81 
cluster that have been corrected for the nebular contami- 
nation so that they give the total stellar spectrum which 
can be compared to the sum of the individual stellar opti- 
cal spectra of the models giving the best UV fits. Fig. ED 
shows that there is again a good agreement between the 
observed and modeled spectra. From the strength of the 
He lines, a "mean" spectral type 07 is found (which is in 
fact similar to the spectral type of the two most luminous 
stars of the cluster, namely star 1 and 2). This indicates 
that despite the difficulties in estimating the stellar prop- 
erties of the individual stars (in particular T c g), our results 
are reliable. Moreover, the He II A4686 line is reasonably 
fitted by the models. As this line is usually filled by wind 
emission, this may be an indication of not too strong mass 
loss rates. 

6.2. Evolutionary status 

We have placed the N81 stars on an HR diagram con- 
structed with Geneva evolutionary tracks without rota- 
tion forZ = 0.004 (Fig. ITI)l . All the stars are compatible 
with an age < 5 Myrs. While star 1 seems to be slightly 
older, the 3 other stars (especially star 3 and 11) lie close 
to the ZAMS and have an age between and 4 Myrs. 



There seems to be an age dispersion of the order 1 or 2 
Myrs in the cluster, which is reasonable. The inclusion 
of rotation is known to move the ZAMS towards lower 
T e g (see Meynet & Maeder 119971) . reducing the age esti- 
mates for stellar populations compared to the non-rotating 
case. However, this effect becomes significant only for rota- 
tional velocities close to the critical velocity. As the SMC- 
N81 stars have V s'mi of only 200-300 km s -1 (which 
corresponds to ~ 1/3 of the break velocity), our results 
are probably not strongly hampered by the use of stellar 
tracks without rotation. 

Based on the equations governing the dynamical evo- 
lution of HII regions fDvson fl978|l with the typical values 
M = 3 x 10" 9 M Q yr" 1 Voo = 1500 km s" 1 , n = 400 
cm -3 (gas density, see Heydari-Malayeri et al. 19.88) and 
a radius of the HII region of 3 pc, one can derive an age of 
~ 1.4 Myrs for the N81 region. This assumes that star 1 
and 2 provide most of the mechanical energy (|Mwoo)- 
Previous age estimates based on the decrement of Hp 
equivalent width by Heydari-Malayeri et al. ( 1988) were of 
1 to 2.5 Myrs assuming a stellar cluster of solar metallic- 
ity. However, for clusters with a small number of ionising 
stars this method is inherently inaccurate due to statisti- 
cal fluctuations. 

If we compare the positions of the SMC-N81 stars with 
the calibration L - T c g for dwarfs of Vacca et al. ^ 1996 ) 
based on studies of Galactic O stars (long dashed line in 
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Fig. 13. Comparison between observed optical spectrum 
of the N81 cluster corrected for nebular lines (solid line) 
and sum of the optical spectra of the models giving the 
best UV fits (dashed line). Each optical spectra from the 
models has been convolved with rotational velocities de- 
rived from the UV spectra. The agreement is very good, 
showing that our determination of the stellar parameters 
of the N81 stars is reliable. 

Fig. 114(1 . we note that except star 1, the other components 
are all underluminous compared to "normal" O dwarfs. 
This result based on the quantitative study of the stellar 
properties confirms the conclusion of paper I in which the 
sub luminosity of most of the stars was already highlighted 
and interpreted as an indication of the youth of the stars. 

In paper I it was argued that the SMC-N81 stars could 
belong to the class of Vz stars. Fig.ll4lshows that although 
young, the stars do not lie perfectly on the ZAMS. This 
may be an indication that either Vz stars are not strictly 
ZAMS stars but more probably "young" stars less evolved 
than "normal" dwarfs, or that our objects are not "true" 
Vz stars (we recall that we do not dispose of optical spec- 
troscopy to firmly establish if the SMC-N81 stars are Vz 
stars or not) . A deeper investigation of the properties of 
this interesting class of objects will be presented in a forth- 
coming paper. 

7. Puzzling wind properties 

7.1. Comparison to previous observations and 
theoretical predictions 

One of the most surprising feature of the spectra of the 
SMC-N81 stars is the shape of the wind lines (mainly 
N v A1238,1242 and C iv AA1548,1551 ): while they show 
no emission, they display a blucshiftcd absorption pro- 



Fig. 14. HR diagram for Z = 0.004 without rotation. 
Different Geneva evolutionary tracks for various masses 
are indicated, together with the ZAMS and isochrones for 
1,2,3,4 and 5 Myrs (data from Lejeune & Schaercr 2001). 
The filled squares give the position of the SMC-N81 stars 
with the typical errors on their position. The calibration 
log L — logTcff of Vacca et al. I|199ti[) for dwarfs is also 
shown by the long dashed line. Note the underluminosity 
of most of the SMC-N81 stars. 



file. This is quite puzzling: how can we produce such an 
absorption without any emission ? Theoretically, in an ex- 
tended atmosphere, the absorption in the P-Cygni profile 
comes from the removal of photons from the observer's 
line of sight, and the emission is due to photons that have 
been scattered isotropically on sight lines not parallel to 
the observer's one. Globally, the amount of "absorbed" 
photons is then equal to the number of "emitted" pho- 
tons. If the atmosphere is only weakly extended (i.e. its 
height is smaller than the stellar radius), then roughly 
half the photons are backscattered towards the star and 
destroyed, so that the emission is reduced compared to 
the absorption (the ratio between them being 1/2). This 
is even more true if there is an underlying photospheric 
absorption which increases the absorption with respect to 
the emission. This points to the direction observed in our 
spectra, but the absence of emission remains a puzzle. 

A possible explanation for the absence of emission 
may be an enhanced backscattering of photons towards 
the photosphere, so that the emission is even more re- 
duced compared to the absorption. But the physical rea- 
son for such a process is not known. Another possibility 
is a strongly non spherically symmetric wind. Rotation 
is known to increase the mass loss near the poles and to 
reduce it at the equator (Bjorkmann & Cassinelli 1993 
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Macdcr & Mcynct 2000). If a star rotates sufficiently fast 
to induce a strong contrast between the polar and equa- 
torial ejection, and if we can observe it pole-on, then we 
can obtain the kind of profile we have: the blueshifted 
absorption comes from the high density polar flow while 
the (absence of) emission corresponds to the low density 
"equatorial" atmosphere. However, as all the SMC-N81 
stars show such profiles and as the probability to see all 
of them pole-on is low, this explanation is not satisfying. 
Moreover, as the projected rotational velocities of the N81 
stars are of the order 200/300 km s^ 1 , it is not likely that 
they are seen pole-on. 

In spite of this curious features, we have been able to 
place constraints on the wind properties of the SMC-N81 
stars, the main results being that they lose mass at an ex- 
tremely low rate. Whether this is typical or not of O dwarfs 
at low metallicity remains to be established. Indeed, this 
is one of the first quantitative determination of mass loss 
rates for such objects and consequently few comparisons 
are possible. To our knowledge, the only previous studies 
of wind parameters O dwarfs in the SMC have been per- 
formed by Puis et al. <|199f>j) . Bouret et al. H2003J) and very 
recently Massey et al. I|2004|) . Puis et al. ifTMB) have only 
derived upper limits on the mass loss rates for four SMC 
dwarfs. This is mainly due to the fact that their method re- 
lies on the the strength of the H a emission whereas the H a 
profile in SMC dwarfs is mostly in absorption. Their upper 
limits are of the order 10 -7 M Q yr _1 which is roughly 10 
times higher than our estimates for the SMC-N81 stars. 
The recent work by Bouret et al. 1|2003[1 include 5 SMC 
O dwarfs of which three (with spectral type between 09.5 
and 06) have mass loss rates between 10~ 10 and 3 x 10~ 9 
Mq yr _1 . Our results agree very well with these estimates 
(we recall that the SMC-N81 stars have mid to late spec- 
tral types). The dwarfs studied by Massey et al. (2004) 
are more luminous than the N81 stars and have therefore 
higher mass loss rates. Fig. 1151 shows mass loss rates as a 
function of luminosity for Galactic and SMC dwarfs from 
different studies. The results of this work and of Bouret 
et al. |QT)03) indicate a strong reduction of the mass loss 
rate of stars with log 4 s - < 5.5. 

How do these results compare to theoretical calcula- 
tions ? The most recent predictions of wind parameters 
as a function of metallicity are the extensive calculations 
by Vink, de Koter & Lamers (2001). Using their cooking 
recipe to estimate the mass loss rate as a function of stel- 
lar parameters and metallicity, we found M of the order 
lO" 7 - 8 M yr _1 for the SMC-N81 stars (the exact values 
are given in TableQJ. The derived mass loss rate are more 
than 30 times lower. Even if we take into account a possi- 
ble error in the ionisation fraction of the CMFGEN models 
(which would lead to an increase of the M determinations 
by a factor of 10) the mass loss rates are still lower than 
predicted (with perhaps the exception of star 3 for which 
we recall that the M determination is uncertain). There is 
thus no doubt that the winds of the SMC-N81 stars show 
an unusual weakness. 
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Fig. 15. Mass loss rate as a function of stellar luminos- 
ity for O stars. Filled (open) symbols are Galactic (LMC, 
small red / SMC, blue) objects. Triangles (squares, cir- 
cles) are for luminosity class V (III, I). Crosses are the 
SMC stars of Bouret et al. I|2003[) . Data are from Puis et 
al. 1(1355)1 . Herrero et al. l|2UUU|> . Crowther et al. lj2UU2j) . 
Repolust et al. (gUU3|), Hillier et al. l|2UU5|l and Massey et 
al. (|2"UU4) . The star symbols are for the SMC-N81 stars. 
Note the low mass loss rates of the SMC objects with 
log -rj- < 5.5 and of 10 Lac (large filled triangle). 



This result is even more striking in term of modified 
wind momentum. This quantity defined by Mvoo^/H is 
predicted to be a power law of the sole stellar luminosity 
(e.g. Kudritzki & Puis. 2000) to give the so-called modified 
wind momentum - luminosity relation (hereafter WLR) 



Rot L- 



(4) 



Fig. ^j] shows the observed relation for Galactic and 
MC stars. For Galactic stars, the correlation is well defined 
and is slightly different for luminosity class I stars and 
LC class III-V stars (as demonstrated by the regression 
curves). For the MC objects however, the relation seem 
to be different from that followed by the Galactic objects, 
at least for the LC V stars. Indeed, our results indicate 
that the SMC-N81 stars (star symbols if Fig. have 
modified wind momenta lower by ~ 2 orders of magnitudes 
compared to Galactic objects of the same luminosity and 
following the LC III-V relation, which simply reflects the 
weakness of the winds 5 . The study of Bouret et al. f 2003 ) 



5 To compute the modified wind momenta of the SMC-N81 
stars, we have adopted Voo = 2.6 v esc since we have not been 
able to derive the true terminal velocities due to the weakness 
of the winds. 
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Fig. 16. Modified wind momentum - luminosity relation. 
Data and symbols are the same as in Fig. The SMC- 
N81 stars (for which we have adopted v oc =2.6 v esc to com- 
pute the modified wind momentum) confirm the trend of 
either a reduced modified wind momentum at low lumi- 
nosity or a steeper slope of this relation at low metallicity. 
The small star symbols show the position of the SMC-N81 
stars if M was systematically underestimated by a factor 
10. The long - dashed (short - dashed) line gives the mean 
relation for luminosity clas V - III (I) stars from Repolust 
et al. H2003JI . The position of 10 Lac is also indicated by 
the large filled triangle. 



indicate a similar trend. Note that even if the mass loss 
rates are underestimated by a factor 10 (see Sect. I5.2.2|l 
there is still a significant difference compared to Galactic 
objects (except in the case of star 3) as shown by the 
small star symbols in Fig. We have also compared 
the derived wind momenta with those expected for B and 
A supergiants (e.g. Kudritzki et al. H999): it turns out 
that the SMC-N81 momenta are weaker. This confirms 
the extreme weakness of the winds of the SMC-N81 stars. 

7.2. Possible origin of the low M /wind momentum 

Several possibilities can be invoked to explain the low 
wind momenta observed for the SMC-N81 stars: first, 
there may be a breakdown of the modified wind momen- 
tum - luminosity relation at low luminosity for dwarfs; 
second, metallicity can affect this relation; third, ion 
runaway may happen in the wind, reducing the mass 
loss rate. In the following, we investigate these different 
hypothesis. Finally we will speculate on the possible 
relation between these "unusual" wind properties and the 



youth of these stars. 

1) Breakdown of the wind momentum - luminosity re- 
lation at low luminosities: 

Let us first consider the case of Galactic stars. The 
study of Puis et al. (1996) showed that in the luminos- 
ity range 5.2 < log L/Lq < 5.5 several of their stars 
seemed to indicate a weaker modified wind momentum 
than expected from the relation followed by higher lumi- 
nosity stars indicating a possible break in the slope of the 
WLR. However, the recent reanalysis of these stars by 
Repolust et al. |2003) do not show this trend any more 
if their upper limits are taken as the true values for M. 
As their mass loss determinations are based on H a fitting 
and as this line becomes insensitive to M below ~ 10~ 8 
M yr _1 (e.g. Herrero et al. I2002|l . one can not rule out 
the possibility of low wind momenta for these low lumi- 
nosity objects. UV analysis of these stars could certainly 
shed more light on this question. 

Concerning low metallicity objects, few LMC stars 
have been studied so far, all of them being furthermore 
of high luminosity. More SMC stars have been stud- 
ied (including the SMC-N81 objects). For objects with 
logL/LQ > 5.5, the wind momenta are lower than for 
Galactic stars, but the WLR seems to have the same slope. 
The low luminosity stars show a clear reduction of the 
wind momenta which imply a steeper slope of the WLR 
in this luminosity range. Nonetheless, most of the high 
L objects are giants or supergiants while low L objects 
are all dwarfs. As the WLR may be different for differ- 
ent luminosity classes, the existence of a real breakdown 
is unclear. 

Moreover, reasons for such a possible break down of the 
wind momentum - luminosity relation are not known. Is 
it linked with the driving mechanism ? A possibility could 
be a change of the ions responsible for the radiative ac- 
celeration with the consequence of the modification of the 
efficiency of the driving, but this would be linked to the 
change of the ionisation in the atmosphere and then more 
to T e ff than to L. Is this break down of the wind momen- 
tum luminosity relation only due to the low luminosity of 
the stars ? There are in fact other objects, central stars 
of planetary nebulae, which are much less luminous than 
any O stars and which seem to follow the mean relation for 
Galactic objects with log > 5.5 (see Kudritzki & Puis 
2000). This renders the behaviour of the low luminosity O 
stars even more puzzling. 

Another possibility is that the current formalism of 
the radiation theory may be erroneous in the case of mas- 
sive stars with weak winds. Indeed, Owocki & Puis ([1999) 
have shown that in such winds, the curvature of the ve- 
locity field near the sonic point could lead to an inward 
directed flux of the diffuse radiation field that can signifi- 
cantly reduce the total radiative acceleration compared to 
the CAK approach (see their Fig. 7). As a consequence, 
the mass loss rates are also reduced compared to the pre- 
dictions based on the CAK formalism for the radiative 



Fabrice Martins et al.: Young massive stars in SMC N81 



19 



acceleration. For stars with high luminosity and/or low 
gravity, the density scale height just above the photo- 
sphere is higher 6 so that line thermalisation near the 
sonic point may suppress the above effect (see Owocki & 
Puls ll995jl . This may explain why giants/supergiants and 
dwarfs with high luminosity are not sensitive to the ef- 
fects of the curved velocity field and show a more classical 
behaviour. This possibility is attractive since it could ex- 
plain why only dwarfs with low luminosities seem to have 
wind properties deviating from the predictions of the CAK 
theory. 

Puis, Springmann & Lennon (2000) have made a de- 
tailed study of the line statistics and its effect on the radia- 
tive driving and have shown that under certain conditions, 
the parameterisation of the radiative acceleration with the 
CAK formalism (g cx (^) Q ) is not valid, and then the pre- 
dictions based on this formalism are erroneous. This would 
happen if both the level density of the ions responsible for 
the driving were much higher and the distribution of the 
oscillator strengths were much steeper. However, they ar- 
gue that the atomic physics of the driving ions should not 
lead to such extreme conditions. 

In summary, the above discussion presents evidence 
that the current predictions of the radiation driven wind 
theory may fail to explain the winds of low luminosity 
O dwarfs due to subtle effects negligible in the winds of 
stars studied so far. 

2) Metallicity effect: 

The wind momentum - luminosity relation may be dif- 
ferent at low metallicity. Such a dependence, though quite 
weak, is in fact predicted by the radiation driven wind the- 
ory. The inverse of the luminosity exponent in equation 0] 
(a) is linked to the line strength distribution function. 
But this function depends on metallicity, a being lower 
at low Z (e.g. Abbott 119821 Puis, Springmann & Lennon 
2000 ) . It turns out that the wind momentum - luminosity 
relation should have a steeper slope at sub solar metal- 
licities. The reduction of the number of lines effectively 
driving the wind also leads to a global shift of the rela- 
tion towards lower modified wind momenta. Our results 
tend to indicate a steeper slope of the WLR, together with 
the Puis et al. (1996) SMC objects and the Bouret et al. 
(2003) results. Nonetheless, the real form of the WLR is 

MwooVfioc (M(l-r))f-i£« (5) 

where T is the ratio of radiative acceleration to grav- 
itational acceleration. For solar metallicity, a ~ 2/3 so 
that the dependence on the effective mass M(l — T) van- 
ishes. Values of a < 2/3 brings back this dependence, 
which should lead to a large scatter of the WLR for low 
Z objects and then to a more difficult calibration of this 
relation. Fig. 1171 shows the modified wind momentum - 
luminosity relation for the SMC stars studied so far: the 

6 The density scale height is proportional to R+/M(l — F) 
where 7?* is the stellar radius and V the ratio of electron scat- 
tering to gravitational acceleration 
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Fig. 17. Modified wind momentum - luminosity rela- 
tion for SMC stars. Data are from Puis et al. 1^996), 
Bouret et al. (2003), Massey et al. (2003) an d tne present 
work. Dwarfs (giants, supergiants) are shown by trian- 
gles (squares, circles). The long - dashed (short - dashed) 
curves are regressions for supergiants (giants/dwarfs) of 
the Galaxy (light curves, from Repolust et al. 2003)) and 
the SMC (bold curves). SMC giants and dwarfs may in- 
dicate a steeper slope of the modified wind momentum - 
relation at low metallicity, but the scatter is important 
and more studies are needed to confirm this trend. 

scatter is indeed important. Indicative regression curves 
are given and reveal a possible steeper slope (a ~ 0.4) for 
giants/dwarfs while for the supergiants a may be close to 
2/3. More studies are needed to confirm this trend of a 
steeper slope at lower metallicities. 

However, two points may impede this explanation. 
First, the most recent hydrodynamical simulations do not 
show any change of the slope of the WLR at metallici- 
ties typical of the Magellanic Clouds (Vink et al. 120011 
Hoffmann et al. I2UU2|> . Furthermore, Kudritzki lj2UU2|> 
has shown that only for extremely low metal content - 
10~ 3 Zq - a difference appears, although his calculations 
were made at high luminosities (logL/L© > 6) were most 
of the driving is thought to be done by H and He lines, 
leading to a weaker metallicity dependence of the wind 
properties. Second, some Galactic stars are known to dis- 
play signatures of weaker than normal winds (Walborn et 
al. H995|l but none has been analysed quantitatively so far. 

We have carried out a detailed study of the wind prop- 
erties of one of them: 10 Lac. This 09V star has often 
been considered as a standard dwarf star for its stellar 
properties. We have run models to fit its IUE and opti- 
cal spectra taking the the stellar parameters from Herrero 
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et al. I|2002|) 7 and the CNO abundances from Villamariz 
et al. ((2)02). A microturbulent velocity increasing from 
5 km s _1 near the photosphere up to 100 km s _1 in the 
outer atmosphere was used. Fig. ^| shows the results of 
our best fit model for which M = 2 x 10~ 9 M Q yr" 1 and 
Voo = 1070 km s" 1 . This value of M is lower than esti- 
mates of Howarth & Prinja (|1989l 10" 6 7 M yr _1 based 
on line profile fitting with a Sobolev code), of Leitherer 
fifflfll 10- 6 83 M Q yr" 1 based on H a ) and lower than the 
prediction of Vink et al. lj27!0T}> - 10" 6 2 M yr" 1 - but is 
in agreement with the more recent upper limit of Herrero 
et al. l|2UU2|) based on H a fits (10~ 8 M Q yr" 1 ). The im- 
portant point is that the wind momentum of this star is 
as low as for the SMC-N81 stars (see Fig.UHl). The metal 
content of 10 Lac being near solar (10 Lac belongs to the 
nearby association Lac OBI, and its CNO abundances are 
between 0.5 and 0.9 the solar values - see Villamariz et al. 
2002|, this indicates that metallicity is not uniquely re- 
sponsible for the weakness of the winds observed in the 
present work. 

The above discussion shows that our understanding of 
the metallicity effects on the wind properties of massive 
stars is still partial and that both observational and 
theoretical efforts are needed to improve this situation. 

3) Decoupling: 

Finally, an alternative answer to the puzzling be- 
haviour of the wind properties of the low luminosity O 
stars is that a decoupling between the absorbing ions 
and the passive plasma may occur. In classical radiatively 
driven winds, ions from metals gain momentum from the 
photons they absorb and redistribute this momentum to 
the passive plasma through frictions. However, if the den- 
sity in low enough, these frictions may become inefficient 
so that the absorbing ions are the only species to be accel- 
erated. Due to their low abundances compared to H and 
He (constituting most of the passive plasma), the mass 
loss is greatly reduced. 

Various studies have been pursued to determine 
the conditions under which such decoupling may take 
place (e.g Babel EMS! HMfl Krticka & Kubat I2UUU|) . 
Springmann & Pauldrach ( 1992) have estimated that 
decoupling occur when radiative acceleration is not 
compensated by frictional braking. We have applied their 
equation 15 to the case of our SMC-N81 stars and we have 
found that decoupling is expected to happen only in the 
outer region of the wind where the outflow velocity has 
almost reached the terminal velocity. In the transsonic 
region, where the mass loss rate is set, decoupling is not 
predicted to take place. More recently, Owocki & Puis 
(2000) have given an estimate of the wind velocity at 
which maximal coupling occur as a function of stellar 
and wind parameters. If this velocity is greater than 
the terminal velocity, then the wind can be described 
by a single component fluid. In the SMC-N81 stars, 

7 T e ff= 36000 K, logy = 3.95, and the He abundance e = 
0.09. 
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Fig. 18. Fit of the UV and optical spectra of 10 Lac (light 
curve: observation; bold curve: model). The parameters 
of the models are: T cS = 36000 K, \ogg = 4.1, M = 
2 x 10~ 9 M© yr" 1 = 1070 km s" 1 . A turbulent velocity 
increasing from 5 km s _1 near the photosphere to 100 km 
s" 1 in the outer wind has been adopted. Abundances are 
from Herrero et al. ( 2002J) for He, Villamariz et al. (2002 ) 
for CNO and have been set to the solar value (Grevesse 
& Sauval lT998|) for the other metals. A rotational velocity 
of 40 km s _1 has been adopted. 

this maximal drift is predicted to happen at velocities 
of the order Voo so that decoupling is not expected to 
happen. Even more recently, Krticka et al. I|2003[) have 
estimated the metallicity below which decoupling should 
occur due to the reduced wind density for a given set 
of stellar parameters. In the case of the SMC-N81 stars, 
this limit is of the order 1/100.Z© which is well below 
the metallicity of the SMC. As a consequence, it seems 
that decoupling cannot explain the weakness of the winds. 

After these considerations the most likely explanation 
for the weak winds of the analysed stars is therefore that a 
breakdown of the modified wind momentum - luminosity 
relation exists at low luminosity for dwarfs (1), possibly in- 
dependently of metallicity. A possible mechanism respon- 
sible for such a behaviour may exist (cf. above) although 
no realistic calculations have been done for parameters 
appropriate to the objects considered here. In any case, 
it is not clear what parameter(s) would determine that 
some O dwarf stars have such weak winds in comparison 
to other O stars of similar luminosity. Possible relevant 
parameters could be a higher gravity, higher mass/light 
ratio, or others (such as the presence of magnetic fields). 
From the current limited sample of objects showing such 
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puzzling wind properties we speculate that the low mass 
loss rate is probably intrinsically related to the youth of 
the stars. It is conceivable that we are beginning to wit- 
ness the "onset" of radiatively driven winds in young, still 
somewhat underluminous O stars shortly after their for- 
mation. Further systematic studies of Vz stars and related 
objects with indications of weak winds will be necessary 
to resolve these issues. 

8. Conclusion 

Based on UV spectral obtained with STIS/HST we have 
analysed the stellar and wind properties of the four main 
exciting stars of the High Excitation Blob SMC-N81 us- 
ing extensive calculations of spherically expanding non- 
LTE line blanketed atmosphere models with the code 
CMFGEN. 

The main results are the following: 

o The stellar properties (L, T e g) indicate that the SMC- 
N81 components are young (~ 0-4 Myrs old) O stars 
which shows, with perhaps the exception of star 1, 
a lower luminosity than "normal" Galactic O dwarfs. 
This, together with the closeness to the ZAMS for star 
3 and 11, confirms the conclusion of paper I that they 
may belong to the Vz class (Walborn & Parker 1992). 

o The UV spectra of the N81 stars show unusually weak 
stellar winds. The upper limits on mass loss rates are of 
the order a few 10 -9 M Q yr" 1 which is low compared 
to 1) Galactic stars of the same luminosity and 2) the 
most recent predictions of M as a function of stellar 
parameters and metallicity. Point 1) could be quali- 
tatively understood due to the reduced metallicity of 
the SMC but point 2) indicates that this reduction is 
higher than expected. 

Although the mass loss rates derived from the UV 
line analysis are potentially affected by uncertainties 
in the modeled ionisation fractions, various tests indi- 
cate that the above conclusions remain qualitatively 
valid. 

o Our objects show modified wind momenta 
{MqVooR 1 / 2 ) which are, for the same luminosity 
L, lower by typically two orders of magnitude com- 
pared to the "normal" O star samples. Similarly low 
wind momenta have also been found by Bouret et al. 
(2003) for 3 SMC stars in NGC 346. 
The modified wind momentum - luminosity relation of 
all the SMC objects could be interpreted as showing 
a break-down at low luminosities or a different slope 
than the Galactic relation. The current sample of SMC 
stars may indeed indicate a steeper slope at least for 
giants and dwarfs, but the scatter is still too large to 
firmly establish this trend. However, the most recent 
hydrodynamical models (Vink et al. 120011 Kudritzki 
2002, Hoffmann et al. I2002f) do not predict such a 
change in the slope between solar and SMC metallici- 
ties. Furthermore we present the first indications that 
some Galactic objects have also low wind momenta 



comparable to the SMC dwarfs. This also tends to ex- 
clude explanations based uniquely on metallicity. 
Possible explanations for a breakdown of the modi- 
fied wind momentum - luminosity relation at low lu- 
minosities are discussed. Ionic decoupling appears un- 
likely according to various estimates. A failure of the 
CAK parameterisation in high density atmospheres, 
discussed by Owocki & Puis l|1999|l . might be invoked 
to explain a lower acceleration in the transsonic region 
where the mass loss rate is set. 

Although the physical mechanism leading to such weak 
winds remains currently unknown, we speculate that 
the low mass loss rate is probably intrinsically related 
to the youth of the stars, possibly testifying of a phase 
of the "onset" of radiatively driven winds in young O 
stars shortly after their formation. 

Further studies of very young massive stars, Vz stars, 
and related objects with indications of weak winds will be 
of great interest to attempt to understand these puzzling 
wind properties and to provide interesting constraints on 
the development of stellar winds in the early phases of 
massive star evolution or possibly even on the final phases 
of their birth. 

Appendix A: Best fits 

Fig. EH El El and El give the best fits achieved for 
star 1, 2, 11 and 3 respectively. For each figure, the model 
parameters are given and the main interstellar lines are in- 
dicated. The interstellar CIV absorption has been added 
to the synthetic spectra as described in Sect.0J Even with 
this correction, the fit of the C IV AA1548,1551 line remain 
poor for star 2, showing the difficulty to produce a signif- 
icant absorption without emission. The normalisation of 
the observed spectra below 1200 A is very uncertain so 
that any comparison with models in this wavelength range 
is irrelevant. 
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